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SUPERSONIC WAVE DRAG OF NONLIFTING DELTA WINGS WITH 
LINEARLY VARYING THICKNESS RATIO 
By Arthur Henderson, Jr. 


SUMMARY 


The supersonic wave drag of a nonlifting, symmetrical, double- 
wedge-profile,, delta wing the thickness ratio of which varies linearly 
in the spanwise direction is calculated by means of linear theory. 

In general it is found that a delta, wing with linearly varying 
thickness ratio can have less wave drag than a constant-thickness-ratio 
delta wing of the same plan form when both wings have either the same 
projected frontal area or the same internal volume. The thickness 
distributions for minimum drag and the corresponding values of the ratio 
of the drag of a linearly varying thickness -ratio wing to a constant- 
thickness-ratio wing are found. 


INTRODUCTION 


In reference 1 , Puckett has found the supersonic wave drag of a 
nonlifting, symmetrical, double -wedge -prof ile, delta wing with constant 
thickness ratio. It is shown that the drag coefficient for the delta 
wing may be reduced below the two-dimensional value only for the case 
in which both the leading edge and the ridge line are subsonic, the 
maximum reduction for a given Mach number and semiapex angle being a 
function of the position of the ridge line. 

In the present paper the additional effect on the drag of varying 
the thickness ratio in the spanwise direction is determined with the 
assumption that the thickness ratio varies linearly in the spanwise 
direction, which means that, when the wing is viewed from behind, the 
line of maximum thickness has a parabolic shape., The source -distribution 
method developed in reference 1 is used to represent the variable- 
thickness-ratio wings. 

It is shown that the deviation of the maximum-thickness line from 
a straight line can be represented by a nondimensional parameter. Charts 
are presented from which, if this parameter is known, the drag of a 
variable -thickness -ratio delta wing can be found. Also shown is the fact 
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that a delta wing with a linearly varying thickness ratio caa give less 
wave drag than a constant -thickness -ratio delta wing when "both wings have 
the same projected frontal area or when both have the same internal 
volume. The value of the nondimensional may 1 mum -t hi r tma RR-rH Kt.ri~hut.inn 
parameter, for which the ratio of the drag of the variable -thickness- 
ratio wing to the drag of the constant -thickness -ratio wing is a minimum , 
is calculated for both criteria. 


SYMBOLS 


Af 

a = P e RL 
B 

c r € LE 

b = P 6 le 


projected frontal area 


P(y/c r ) 

b 



/ 


CD 

c r 

c(y) 

D 


Ap 

€ LE 


e RL 

Fi,F 2 ,F 3 

Gi,G 2 ,G 3 


drag coefficient, D/qS 
root chord 

local chord 
drag 

local static pressure minus free-stream static pressure 
tangent of leading-edge semiapex angle 
tangent of ridge -line semiapex angle 

functions defined in appendix B 


X 


slope of airfoil surface in streamwise direction 
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M 

m 

m 

P]_ and p-^ 

P2 and P2 
Ql,Ql and Ip; ip 
Q 2 j Q 2 and q 2 ,q 2 


free -stream Mach number 

v 

dimensional thickness -distribution parameter 

nondimens ional thickness-distribution parameter, 
mc r € LE mc r 8 
t “ Pt 

dummy variables for x and y, respectively 

pres sure -distribution function and 80 p/ 8 x associated 
with basic source distribution when b < 1 

pressure-distribution function and 802 /^ associated 
with superimposed source distribution when a < 1 

pressure -distribution functions and 80 p/ 8 x associated 
with basic source distribution when b > 1 

pres sure -distribution functions and 802/8x associated 
with superimposed source distribution when a > 1 


q 



P 

S 

*r 

t(y) 

T 


dynamic pressure. 



free -stream density 
plan-form area 
root thickness 
local may i mum thickness 
thickness ratio at root, t r /c r 


V 


free-stream velocity 
internal volume 


Vi 

0 


perturbation velocity potential 



h 
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v perturbation velocity in z-direction 

x,y, z Cartesian coordinates (right-handed system), x in 

direction of free -stream velocity 

Subscripts: 


R 

1 

2 


right half of wing 
basic distribution 
superimposed distribution 
Primes denote constant-thickness-ratio wing. 


ANALYSIS 


The linearized partial differential equation for the perturbation 
potential in steady supersonic flow is 


0 


2 8 % 

" «« ■ w l — 

8x^ 3y^ bvfi 


0 


( 1 ) 


When considering a thin symmetrical airfoil at zero angle of attack, 
for which small disturbances axe assumed, the boundary conditions on 
the surface of the wing may be satisfied, to the first order, in the 
plane of the wing. Hence 



( 2 ) 


In reference 1, it is shown that a solution of equation (l) which satis- 
fies condition (2) is that for a source distribution 


_ _i ff w(g,q)dg dq 

\/(x - S ) 2 - P 2 (y - n) 2 ' 


0U,y) 


( 3 ) 
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where w(i,ri) = X(I,tj)V is proportional to the source strength per 
unit area and X is the slope of the airfoil section measured in the 
streamwise direction. For the wing shown in figure 1, the slope of 
the surface ahead of the ridge line is equal to the ratio of one -half 
the local maximum thickness at the ridge line to the part of the local 
chord ahead of the ridge line; the slope of the rear surface is the 
negative ratio of one -half of the local max imum thickness at the ridge 
line to the remainder of the local chord. In the present paper the 
thickness ratio is assumed to vary linearly with y, or 


t kl = 

c-(yj 


+ 2my 


w 


where t r /c r is the root thickness ratio, and m is a parameter which 
dete rmine s the shape of the thickness curve in the spanwise direction. 
Therefore, over the forward part of the wing 


v(x,y) 


MzlZ 

2(1 - r)c(y) 


V / t r 
(1 - r)\2c r 


+ my 


(5a) 


and over the rear part 


w(x,y) 


- t(y)v 

2rc(y) 



(5b) 


where r = and represents the ridge-line position in percent chord 

measured from the trailing edge. 

In the following analysis c$/cSx is evaluated from equation (3) by 
integrating over the right half of the upper surface of the. wing. Inas- 
much as the wing shape, and therefore the pressure distribution, is sym- 
metrical with respect to the x-axis, 80(x,y)/8x for the entire upper 
surface is . 
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dfiKx.y) 


d0(x,-y) 

8x 

8k- 

L Jp 

. 8x 


( 6 ) 


The potential of any point (x,y) is found by integrating over the area 
of the wing lying within the forward-drawn Mach lines from the point (x,y) 
in order to obtain the effect on the point of al 1 the sources that can 
influence it. 

The wing of the present paper is represented by a distribution of 
sources in the plane z = 0 . This distribution is uniquely determined 
by the boundary conditions on the wing. Since equation ( 3 ) is a solu- 
tion to a linear partial differential equation, this distribution may 
be represented by a superposition of source distributions in any manner 
desired as long as the net result is that prescribed by the boundary 
conditions. The drag on the wing is then found by computing the drag 
caused by the pressure field of each source distribution and algebrai- 
cally adding the results. 

In the present analysis two source distributions are used: one, 

the basic distribution' which extends over the entire plan form nnrl 
follows the law governing the source strengths representing the forward 
slopes of the wing, and the other, the superimposed distribution which 
lies in the region bounded by the ridge lines and the trailing edge. 

The sources in this region have the strengths prescribed by the boundary 
conditions for this region plus the negative value of the strengths of 
the first distribution. Thus, when superimposed, the net soupce distri- 
bution for the whole wing is that prescribed by the boundary conditions. 

From figure 2 (a), the basic source distribution used to represent 
a delta wing with subsonic leading edge ( cari expressed as 


-- 


«(1 - r) 



+ mqjdt dq 


x-p(y-q) 

h/sle \/(x - 1)2 _ 32( y _ ^2 
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and the superimposed source distribution when the ridge line is sub- 
sonic (P6 rl < 1), as 



By the use of equations (6) and (7)^ a solution for the basic distribu- 
tion is obtained: - 



p2jt(l - r) (1 - b2) 

= Pi(m,r,x,u,b) (9) 
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where u = h = Pst/e, and a = pepp. Equation (9) is obtained in 

appendix A. S imi larly, equations (6) and (8) produce the solution for 
the superimposed distribution: 


8x 


= p 2 = - i p 1 (m,r,x>u,a) 


( 10 ) 


The expression for the potential of the basic source distribution 
used to represent a delta wing the leading edge of which is supersonic 
(0epE > 1) is similar to equation (7); however, since in this case the 
Mach lines from the apex lie on the wing, there are two regions of inte- 
gration (fig. 2(b)). One is the region interior to the Mach lines where 
the right and left halves of the wing are interdependent and 01 ^ is 

exactly equation (7). The integral is evaluated differently, however, 
since pSr.E > 1. The other is the region that lies between the Mach 
lines and the leading edge where the flow is two-dimensional. Here the 
right and left regions are independent of each other, and the expression 
for the potential on one side also gives the potential for the other. 

For the region of the wing exterior to ,the Mach lines. 


c-P(y-Tl) 




V 

It(l - r) 


I q 

1- 0 6 LE u e LE 


2c, 


+ mq I d £ dq 


V(* - D 2 - P 2 (y - 11) 2 


v 

«(1 - r) 


fLE( x+ Py) n x+ P(y-^) 


,2c. 


+ mq d| dq 


q/e 


LE 


V(x - I ) 2 - p 2 (y - q ) 2 


( 11 ) 


which, with equation ( 7 ), yields for a supersonic-leading-edge delta 
wing, when O^uSl 
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2 X 


. a \ ill 

^ pjt(l - r)v/b2 - 1 Vi - u2 


2V / — £\b 


2Vmb 2 x 


tan 


-in/ 


b^ - 1 


p 2 «(l - r)(b 2 - l) 3 ^ 2 \/l - u 2 


2Vmb^ux 


.-1 u ^] 


P 2 jt( 1 - r)(b2 - 1) 


, tan-1 vV Z Lz l : 
3/ 2 -h\/l - u 2 


2Vmb 2 x 


p2«(l - r)(b 2 - 1) 


\/l - u £ 


= q-L(m,r,x,u,b) 


(12a) 


and, when 1 = u % b 


8x 


•(i> 


Vmb^ux 


P(1 - r) \/b 2 - 1 p 2 (l - r)(b 2 - l) 3 ^ 2 P 2 (l - r)(b2 - 1) 


Vmb^x 


J/2 


= q. 1 (i« l ^,x,u,b) 


(12b) 


The corresponding functions for the superimposed source distribu- 
tion when the ridge line is supersonic are for 0 ^ u ^ 1 


802 p 

&T = q 2 = ~ r 


(13a) 
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and, for 1 u ^ a. 


2 _ _ 1 _ , 

“ q 2 “ " ? qpU^x^a) (l 3 b) 


The parameter m in the preceding equations has the dimensions of 

(length) -1 . It will he useful to examine the curves of thickness plotted 
^•Sainst spanwise distance that can be produced by a variation of m 
( 11 S* 3) with the purpose of fo rmin g a new dimensionless parameter which 
includes m. From equation (4), a nondimens ional form of the thickness 
curve is 


t (y) = 

t r 



and since, from figure 1, 


c(y) 



(14) 



(15a) 


where 


0 ^ 


p (r/ c r) < 


The nondimens ional thickness-distribution parameter is, therefore, taken 
to be mc r b/pT which is hereinafter denoted by m. Equation (15a) may 
be rewritten as follows: 
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+ (2E - 1)B - 2mB 2 ‘ (158) 

_ 3(y/c r ) 

8 

Equation (158) is a quadratic in B and represents a paraBola. If the 
type of curve for equation (158) is known., one absolute limit on m can 
8e esta81ished. This limit can. he determined since, at the wing tip the 
slope of the thickness curve cannot he less than zero. For this case, 

m = -i. When the slope of the thickness curve is zero at the root, 

m = however, this value is not necessarily an upper limit. Curves of 

thickness distribution for different values of m are shown in figure 3- 


where 


*(y) _ 


= 1 


B 


PRESSURE DISTRIBUTION 


To the first order, the pressure on the surface of the wing is 
given by 


Ap _ 2/^A 

q VySyz=0 


(16) 


Introducing the parameter m into equations (9), (10), (12a), (12b), 
(13a), and (13b) and using equation (l6) gives the pressure equations in 
the following forms: 


For the basic suhsonic-leading-edge delta wing 


= 


2 rh 


q Ml - r)\/l - h 2 


tanh 


-1 ~ *> 2 
\/l - u 2 


P*(l " r)(l - 82)1^ _ b 2 


tanh' 


-1 \A ~8 2 

Vl - u 2 


hu 


_1 uv/l - h^ ./ r 
tanh y - VI - U c 

^1 - b^ 8 Vl - u^ 


= P ( m > r >Jp u > 1, )' = P 1 


(0^u<8) (17) 
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For the superimposed subsonic-leading-edge delta wing 

^ = _i p(5,r,^-,u,a) = -| P 2 (0 ^ u < a) (l8) 

For the basic supersonic -leading-edge delta wing 



= Q(m,r,^pu/b^ = Q x (l ^ u < b) (19*>) 

For the superimposed supersonic -leading-edge delta wing 


_ 

q 


HIM 

II 

(0 ^ u ^ 1) 

(20a) 

- 

q 


= -H 

(l ^ u ^ a) 

(20b) 
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The pressure field associated with the basic source distribution 
is shown in figure 4 for both the subsonic- and supersonic -leading -edge 

cases for particular wings. For these wings m = which prescribes 

that the maximum-thickness line, which is at the trailing edge, have 
maximum concavity (see fig. 3)* The effect of having m < 0 is that, 
as x/cp increases, the pressures are progressively decreased in a 
spanwise direction below the values for a conical pressure field. If 
m > 0, the pressures would be progressively increased in a spanwise 
direction above the conical pressures as x/c r increases. 


DRAG 


For the drag of the delta wing herein considered, the following 
three distinct cases must be treated: 

Case I - Supersonic leading edge and supersonic ridge line (b > 1, a >1) 
Case II - Subsonic leading edge and supersonic ridge line (b<l, a>l) 
Case III - Subsonic leading edge and subsonic ridge line (b<l, a<l) 

In computing the drag by the method of superposition, it must be 
kept in mind that the basic source distribution creates a pressure field 
the influence of which is felt over the whole surface of the wing; 
whereas the pressure field from the superimposed source distribution is 
felt only over that part of the wing which lies, between the ridge lines 
and the trailing edge, except for case III where the pressure field 
"spills over" the ridge lines and exerts its influence on the forward 
slopes up to the Mach lines from the ridge -line apex. 

Since the wings are symmetrical about both the x,y- and x,z-planes, 
the drag is computed for the right half of the upper surface of the wing 
and the result is then multiplied by four. 

Case I .- The drag for case I (fig. 5(a)) is represented by 


-= r^xd s (21) 

3 Js l 


which, as a result of the superposition of source distributions, becomes 
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t - i ^ ^ + i + B 2 ^ “ - 

i r 02^2 dS + r 

J B-|_ Ja „ 


dS + 


'*2 


f Ql^2 ^ ” r f ^2^2 

j b 3 j b 2 +b 3 


dS 


( 22 ) 


The l imi ts of integration on the second and fifth integrals make 
the operation rather cumbersome. This difficulty may be obviated by 
rewriting equation (22) in the following form, which includes the sub- 


stitution X 2 = ~ 


1 - r 


U-J f QiXj. dS - r as + 

^ Ja 1 Ja 1 +b 1 +b 2 

— / Qg^l dS + i dS - 

r 2 J B-l J A2 


1 - r 



> 

r J Ag+B^ 

For 

this case 

b > 1, a > 1. 


The element of area is < 

are 

x and u 

where u = 

K x 

of 

x and y 

is taken with ] 


Q-|_X2. dS 


1 - r 


f. 


2+ b 3 


dS 


(23) 


becomes dS = — dx du. Perfo rming the operations indicated in equa- 

ls 

tion ( 23)7 reduc ing the result to coefficient form, and making the sub- 
stitution a = — yields 4 • 
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= F x (r,b) + SF 2 (r,b) + 5 2 F 3 (r > b) 


(b > 1 > r) (2k) 


Equations for F]_, F 2j and F^ are given in appendix B and are 
plotted in figure 6. 

Case II .- The drag for case II (fig. -5(l>)) is given by 

= f ds + r P^ dS - ~ I 02^2 dS - J dS 

^ J A ^'B 1 +B 2 ^Bi ^B~ 


r / p iM. «s 

r J A 


- — / P-, X-, ds + 1 — - r / QpX-, dS + 

r 4 + Bi + B 2 r2 J B i 1 


rl 


2 / ^2 X 1 dS 
* J B2 


(25) 


In evaluating equation (25) it must be remembered that b < 1, a > 1. 
The final value may be expressed in coefficient form by 




= G-^(r,b) + mGg^b) + m^G^r^b) (1 > b > r) (26) 


Equations for G-j^ G 2 , and G^ are given in appendix B and plotted 
in figure 7 with the H-functions which are presented subsequently. 


4q 


Case III .- The drag for case III (fig. 5(c)) is expressed as 


dS 


/ Pi^i dS + j 

r P^ ds - i 1 

r P 2 X 2 as - i / 

■^a 1 +a 2 J 

B J 

B 


“'I 


(1 - r) 


a 1 +a 2 


PiM. dS r 


/ 

vJ A 1 +A 2 +B 




(1 - r) 
r 2 


'B 


i f P 2 \i dS 
r j a 2 


?2^1 — 
(27) 
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Here b < 1, a < 1, and the result is 

— 2 ^. = Hn (r,b) + mHo(r,b) + m^Ho(r^b) (1 > r > b) ' ( 28 ) 

T 2 


Equations for Hj_, and are given in appendix B and plotted 

in figure 7 - 


, Aether a variable-thickness-ratio delta wing is better from the 
wave -drag standpoint than a constant -thickness -ratio delta wing depends 
upon the criterion chosen. Equations (2h) , (2 6 ), and (28) give the drag 
for a const ant -thickness -ratio delta wing when m = Oj the results are 
identical to those of reference 1. Inasmuch as F 2 > F 3 , G 2 > G 3 , 
and Hp > H 3 always , it is obvious that, if the variable- and constant - 
thickness-ratio wings are compared on the basis of identical plan form, 
r,b, and thickness ratio at the root, the concave thickness-distribution 


curve 


(-I s 


m 


< 0, see fig. 3) will give lower drag. This criterion. 


however, is a poor one because the internal volume and the projected 
frontal area are reduced. Since the same projected frontal area appears 
to provide a fair basis for comparison and internal volume is important 
in practical design considerations, these two criteria are developed for 
the purpose of comparing the drag of a constant- with a variable-thickness- 
ratio delta wing. 


PROJECTED -FRONTAL -AREA CRITERION 


If a variable -thickness -ratio delta wing and a constant -thickness- 
ratio delta wing have identical plan form, r, b, and projected frontal 
area, then 



c r e LE 


t(y)dy 


and from equation ( 15 b) 


t(y) 


= t T 



2m - 1 


2m 



e LE 2 



c r e LE 


y - 
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Therefore, 


Af = Sr^l + | m) (29) 

where A^ is the projected frontal area, S is the plan-form area, 
and t is the thickness ratio at the root. 

For a constant-thickness -ratio delta wing 

Af * = S't' (30) 


From equations (29) and (30), the stipulation that Af = Af' and 
S = S' gives the relationship between t' and t as 

T' = r(l +|l) (31) 

From equation (2k) 

CpP _P 

— = F x + mF 2 + mF 3 


for the variable-thickness-ratio delta wing, and 


VP _ F 

^' Fl 


(32) 


for the constant-thickness-ratio delta wing. Forming the ratio of equa- 
tion (2k) to equation (32)' and substituting equation (31) gives: 


Cp F-j_ + 5F 2 + m% 3 



( 33 ) 
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The ratio Cd/Cd' can be minimized with respect to m, and the value 
of m, which gives minimum drag for the conditions stipulated, is 


m 


bF ± - 3F 2 

6F3 - 2F 2 


(3b) 


Equation ( 34) applies to the G- and H-f unctions as well- as to the 
F-functions where the corresponding G- and H-functions are substituted 
for the F -functions. 

The nondime nsional thickness -distribution parameter m for minimum 
drag, based on frontal -area considerations, is presented in figure 8. 
Figure 9 presents Cd/Cd' as calculated from equation (33) using the 
values of m obtained from figure 8. When b > 1 > r, it can be seen 
from figure 8 that the wing for minimum drag is essentially that with 
constant thickness ratio (m *** 0), except when b = 1.001, and even for 
this case figure 9 shows a maximum drag reduction of less than 2 percent. 

Men 1 > b > r, figure 8 indicates that the wing for minimum drag 
has a concave thickness distribution (m < 0). The corresponding curves 
of figure 9 indicate a maximum drag reduction of only about 4 percent 
for this case. The main conclusion to be drawn then is that, with the 
projected frontal area held constant, the thickness can be shifted from 
the tip to the root with no penalty in wave drag. 

When 1 > r > b, the optimum thickness distribution is sometimes 
concave (m < 0, for which the same conclusion can be drawn as for 
1 >b > r) and sometimes convex (m >0). Some of the indicated shapes 

for m > 0 are highly unrealistic ^that is, m » 7^ which indicates 

that, within the realm of physically practical shapes, there is no 
minimum m. This condition, however, does not mean that a drag reduc- 
tion is not possible. Figure 10 presents the curve of Cd/Cd' plotted 
against m for a case when m, as determined from equation (34), gives 
an unrealistic shape. It can be seen that, although m for minimum 
drag as calculated from equation (34) is 3-68, values of m from l/2 
to 1 give a drag reduction of 20 to 25 percent. In figure 9 the dashed 
lines show values of Cd/Cd' when m represents a wing that is unreal- 
istic, where an unrealistic wing has arbitrarily been 1 chosen to be one 
for which m > 1. (See fig. 3-) 
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INTERNAL -VOLUME CRITERION 


In order to compare the variable- and constant-thickness-ratio 
delta wings on the basis of internal volume, both wings are specified 
to have identical plan form, r, b, and internal volume. The internal 
volume of the wing of figure 1 is 


Vi 


-s 


z dx dy 


(35) 


For the right half of the wing, from the leading edge to the ridge line 


z 



(36) 


and from the ridge line to the trailing edge 



(37) 


Using equations (35), (36), and (37) gives, for a variable-thickness- 
ratio wing. 


V 


i 



(38) 


and for a constant-thickness-ratio wing 


v - 


T'C, 


(39) 


When equations ( 38 ) and (39) are equated and the conditions that 
Vi = Vj[', c r = c r ', and S = S' are fulfilled. 
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(40) 


Forming the ratio of 
gives 




as before for the present criterion. 


Cjj + mF2 + m^F^ 



(41) 


which, when mi n imized with respect to m, gives 


m = 2 



(42) 


for the value of m for minimum drag under the conditions stipulated. 

Values of m for minimum drag based on internal -volume considera- 
tions are presented in figure 11. Figure 12 presents Cp/Cp 1 as calcu- 
lated from equation (4l) by the use of the values of m from figure 11. 
When b > 1 > r, figure 11(b) shows that the value of m for minimum 
drag is approx ima tely constant at about m = -O.V?j that is, the thickness 
distribution for minimum drag is very close to that of maximum concavity. 
The corresponding curve of figure 12 gives a maximum drag reduction of 
between 8 and 10 percent. 

When 1 > b > r, figure 11(a) indicates a value of m which is 
fictitious in that m < - ^ which would prescribe a wing where part of 

the upper surface is below the bottom surface and conversely. Therefore ,- 
there is no minimum drag for any real wings in this range. Figure 13 is 
a plot of Cp/Cp ' against m for a particular configuration in this 
range. Although the optimum m as determined from equation (42) is 

fictitious, it can be seen that any value of m between m = and 

2 

m = 0 will give a drag reduction, the reduction in drag at m = -i 

being about 20 percent. In figure 12, the values of Cp/Cp' which were 
obtained by using fictitious values of m are represented by short 
dashed lines. 
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When 1 > r > b, figure 11 ^hows that the m for minimum drag 
varies from negative to positive with 'increasing r for each value of b. 

Included in this range are both fictitious ^m < and unrealistic 

^m » 7 ^ values of m calculated from equation (42).' The corresponding 
values of Cp/Cp 1 are shown in figure 12 where the short dashed lines 
were obtained by the use of fictitious values of m ^m < - i j and the 

long dashed lines indicate the use of unrealistic values of m (m > 1 ). 

It should be emphasized, however, that substantial drag decreases can be 
realized in the regions of the dashed curves for practical wings, as 
can be seen from figures 10 and 13 * which indicate the trend of all curves 
of Cp/Cp 1 plotted against m when the m for minimum Cp/Cp ' is 
either positive or negative. 


ILLUSTRATIVE EXAMPLE 


If, for a given constant-thickness-ratio delta wing, it is desired 
to increase the thickness of the root and, for structural reasons, the 
thickness at the tip is not to be decreased, two methods are suggested. 
Either a variable-thickness-ratio wing with prescribed root and tip thick- 
ness ratios or a new constant-thickness-ratio wing with prescribed root 
thickness ratio could be formed. 

For the former case, if t ' is the thickness ratio of the original 
constant -thickness -ratio wing, the slope of the maximum thickness line 
with respect to y is -t '/26t.f. If the variable -thickness-ratio wing 
has this slope at the tip, equation ( 15 ) gives 


a 



y= c r e IE 



+ 2m) 



or 


t' = t( 1 +-2m) 


(43) 


If the root thickness of the original wing is increased by n times. 
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T = XLT 1 (44) 

The two equations (43) and (44) -yield the necessary thickness distribu- 
tion 



(45) 


As an example, consider the wings for which b = 0.8 and r = 0 . 5 . 
the root thickness of the original constant-thickness-ratio wing is 


increased by l/ 2 , n 


3 

2 


and 



When the ratio 


If 


Cjq G^_ + mG 2 + m^G^ 

c d' (1 + ^) 2 G 1 


is formed, the drag of the variable-thickness-ratio wing is found to be 
1.77 times that of the original constant-thickness-ratio wing. 

If the increased root thickness is attained by forming a new 
constant-thickness-ratio wing, the drag is 2.25 times that of the 
original wing. Here, of course, both the projected frontal area and 
the internal volume are larger than those for the variable-thickness- 
ratio wing. The drag of a constant-thickness-ratio wing which has the 
same projected frontal area as the previously mentioned variable-thickness- 
ratio wing is I .78 times the drag of the original wing; whereas, the 
drag of a constant-thickness -ratio wing which has the same internal 
volume as the previously mentioned variable -thickness-ratio wing is I .89 
times that of the original wing. 


CONCLUDING REMARKS 


The variable -thickness -ratio delta wing has been compared with the 
constant-thickness-ratio delta wing under the conditions that they have 
identical projected frontal area and identical internal volume . For 
these conditions the optimum value of the thickness-distribution parameter 
has been determined. The wing shape given by this value, although no t 
always practical or real, nonetheless gives an indication of the types 
of maximum-thickness-line distribution; that is, whether it should be 
convex, straight, or concave. From the analysis the following conclusions 
are drawn: 
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1. On the "basis of fixed frontal area, 

(a) When the leading edge and ridge line are both supersonic, 
the constant -thickness -ratio delta wing has essentially optimum 
wave-drag characteristics for all combinations of lea ding -edge 
semiapex angle and position of the ridge line. 

(b) When the leading edge is subsonic and the ridge line is 
supersonic, the concave maximum -thickness distribution has the 
optimum wave-drag characteristics which, compared with a constant- 
thickness -ratio wing, gives greater thickness at the root and less 
thickness at the tip. It should be noted, however, that the drag 
reduction which can be realized in this range is small. 

(c) When the leading edge and ridge line are both subsonic, 
as the ridge line moves forward from the sonic condition to the 
leading edge, the maximum-thickness distribution for optimum wave 
drag goes from concave to straight to convex. For the latter case, 
the drag decrease indicated is large and absurd values are approached 
as the ridge line approaches the leading edge. 

2. On the basis of internal volume considerations, ■ 

(a) When the leading edge and ridge line are both supersonic, 
the concave maximum-thickness distribution has optimum wave-drag 
characteristics, the drag reduction being about 8 to 10 percent for 
all combinations of leading-edge semiapex angle and position of the 
ridge line. 

(b) When the leading edge is subsonic and the ridge line is 
supersonic, the concave thickness distribution has the best wave- 
drag characteristics, with drag reductions of as much as 20 percent 
being predicted. 

(c) When the leading edge and ridge line are both subsonic, 
as the ridge line moves forward from the sonic condition to the 
leading edge, the maximum -thickness line goes from concave to 
straight to convex $ whereas the corresponding drag reductions range 
from about 20 percent to zero to absurd values. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., October 1, 1952. 
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APPENDIX A 

EVALUATION OF /pe^ < 1) 


Equation ( 7 ) is 


(0l). = 


,y ~x-p(y-T]) / t T 


+ mT] di dT > 


jt(l - r ) J q Jrj/e Lg \/(x - £) 2 - P 2 (y - t ]) 2 
6 LE( X+ Py) px+PCy-n) / t \ . 

te - + mT >) dl dTi 


rt(l - r) 


1+ P 6 le 

U *l/ € LE \/(x - £) 2 - p2(y _ q)2 


Let 


Then 


x - I 


P(y - n) 


= cosh z 


(01)* = - 


V 


R "it(l - r) 




, g LE(x+py) 

1+ P e LE / t. 


x- JQL_ 
e LE 


P(y-R) 

dz - 

cosh - - 1 ' 

x-^L. 

LE 

P(y-R) 


ic(l - r) 


v 2 c 7 


+ mq] cLt] 


dz 


'O 


-^f te * -) 


cosh' 


-1 


x — — 3 — 
g LE 


P(y - t) 


gLgU+Py) 

V_ I l+S*tE _ (Jr + \ cosh -l 
*(l " r ).L \ 2c r / 


dq 


_ _D_ 

e LE 


P(y - r) 


dq 
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Differentiating 


(0i)r with respect to x 


gives 



Ve LE 
ir(l - r) 


=LE 


(x+py) 


1+ P e LE 



By the use of equations (l6o) and ( 169 ) in reference 2, the following 
equation is obtained: 



V 

w 

> 

I 

II 

** 

W y 

R «(1 - r)^l - P^le 2 

2c r 


m (P 2 y e LE 2 - x€ Le) 

lOff 

p(y - xe LE) 

1 _ n2,. 2 

x P € le 


]/{x 2 - p S 2 ) (l - P 2e LE 2 ) + P 2 ^ - x 


m \/* g *LE 2 - 

\A - e 2 *™ 2 


Let b = pe-^ and u = p^q then 



Ve LE 

*r 

) 

R «(1 - r) \/l - b 2 

2c r 


m ( bux eEE " 


l°ge 

b - u 




1 - 



1 - bu - \/(l - u 2 ) (l - b 2 ) 


^LE \A ~ 

\/l - b 2 _ 
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Using equation (6) gives: 
80! Ve 


LE 


8x 


Jt(l - r)\Jl - b 2 


2c, 


l°ge 


b 2 - u 2 


(\/l - u2 - Vl - b2) 


HDC6 


LE 


1 - b 
mbuxe 


logs 


b^ - u c 


(\/i - u2 _ \/rT^ 


2mx6 T - F ,yi - u 2 
\/l - b2 


LE 


1 - b2 


i°ge 


(b - u) 

1 + bu - \ j\ 

[l " ^ 2 ) 

(l 

(b + u) 

1 - bu -\f' 

(l - u2) 

(l - t2)_ 






Finally, after some algebraic manipulations are made and the identity 

(x 2 < 1) 


logg 1 + x = 2 tanh - ^- x 


1 - x 

is used, this equation is obtained in the following form: 

/ t. 


80i 

8x 


2T ls3 b 


tanh' 


-1 £ 


p*(l - r)]/l - b 2 \/l - u 2 ' 


2Vmb 2 x 


3/2 

p2jt(l - r) (1 - b2) 

2Vmb3ux 

3/2 

p2n(l - r) (1 - b2) 


tanh 


tanh' 


-1 \A - ~b 2 _| 

l/l - u 2 

_1 u\/l - b 2 
b\^L - u 2 


2Vmb 2 ? 


Vl - u 2 


p2jt(l _ r)(l - b 2 ) 
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Fl(r/b) 


F 2 (r,b) 


F 3 (r,b) 


APPENDIX B 

PRESENTATION OF F-, G-, AND H-FUNCTIONS 
F-Functions (b > 1 > r) 


(l - r 2 )/ 


-1 1 

COS q- + 


2b , 


b 2 - 1 


(l - r 2 )/ 


tan' 


-1 7b 2 - r 2 


b 2 - r 2 


b - r 


3* 


b(l + r) b(l - r) 


r(l - r 2)i/t2 - r 2 2(1 - r 2 )(b 2 - r 2 )' 


tan 


■1 Vb^Tr 2 


b - r 


b(l - r) _1 1 . 

- ~ ^ - 3/ 2 COS b + 


2(1 - r 2 )(b 2 - l) 


b 2 - r 


2b (b 2 - l)(b 2 - r 2 ) 2r ^ _ r y b 2 _ r 2 


-1 r 

cos M 7- 
b 


k 

n 


b(l - r)(l + Xr 2 ) _ gb(l - r) 2 fl ± 3r 2 ) 


i(l - r 2 ) 2 (b 2 - r 2 ) 3/2 3r(l - r 2 )Vb 2 - r£ 


tan' 




b - r 



b(l - r) 

I 3 + r 2 ] 


6 

[l - r 2 ) 2 (b 2 - l] 

3/2 


b(l - r) 2 

(3 * 4 

3(l - r 2 ) 

3 /b 2 - 1 


1 1 

COS"- 1 - — + 


b 2 - r 


h 6b (b 2 - l) (b 2 - r 2 ) 


2b(l - r) 2 (b 2 + r) jtbfeb 2 - 3r 2 ) 

n ' , 




3 (l - r 2 ) (b 2 - l)(b 2 - r 2 ) 6r(l - r)(b 2 - r 2 ) 
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Gi(r,b) 


G 2 (r,b) 


G-Functions (l > b > r) 


2 

- -< 
n 


(l - r 2 )\/l - b 2 


loge b + loge 


\Il - b 2 


2b 


tan 


-1 \/b 2 - r 2 


1 _i 

+ — 7 - r COS b 


r(l - r 2 )\/b 2 - r 2 1 - r + l/l - b2 ^ 


8 

3 * 


(l - r2) (l - t2) 3 / 2 


log b + 
e 


b(l - r) 1q b 

a(L - -* 2 ) 3/2 l-^i + 


b(l - ri 


(l - r 2 )(b 2 . r2) 3 / 2 
b(l + r) 


r(l - r 2 )\/b 2 - r 2 


tan 


-i 'fcF - f 


2 „2 


1 - r 




b 2 - r 


(1 + r)(b 2 - r) 7 ^ + 

2b (l - r 2 ) (b 2 - r 2 )v^T^ _ 2b ^ ' b 2 ) 


2 r(l X - r) a 03 " 1 b + 


2r(l - r)'/b2 - r2 


cos 


-1 £ 
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G 3 (r,b) 




k 


4b (1 - r) 2 i b3(l - r) 2 



3(1 - r 2 ) Vl - b 2 3(l - r a ) 2 (l - b 2 ) 3/2 


b(l - r) 2 (3 + r 2 ) + 

3(1 - r 2 ) Vl - b 2 


b(l - r) (3 + r 2 ) 

1 nr ^ 1 

6(l - r 2 ) 2 (l - b 2 ) 3 ^ 

Xo 8e + 

1 - \/l - b 2 

2b(l - r) 2 (3 + r 2 ) + 


3(1 - r 2 )Vb 2 - r 2 


b(l - r) 2 (l + r 2 ) 

tan -1 ^ ~ r2 1 

3(1 - A) a (b2 . r2) 3/2 

J 1 - r + VI - b 2 

(1 - r) [a(l - b 2 ) - (t 

)2 ~ r )j 

6b (l - r 2 )Vl - b 2 

r(l - r) (l - b* 4 ) 

b 2 - r 


6b(l - r 2 ) 2 (b 2 - r 2 ) 6»(* 2 - r 2 ) (l - ^ 

2b(l - r) 2 (b 2 + r) + bfeb 2 - 3r 2 ) 

3(1 - r 2 ) 2 (b 2 - r2)(l - b 2 ) 6r(l - r)(b 2 - r 2 ) 3 ^ 


o'l'i 
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Hp^b) 


H2(r,b) 


2 

it 


(l - r2)v/l - b2|_ 


H-Functions (l > r > b) 


logg b + logg 


b 

1 - \/l - b 2 __ 


log. 


(1 - r)l 

[r -f V'r 2 ' 

1 

o' 

rol 

i 

r - b2 + \J\ 

(r2 - b2) 

(x- 

l 2 ) 


1 . _i . 1 , 4 b 

sin b + sin — - 


r(l - r) 


r(l - r) 


b , b 

p--— loge - + 

•(l - r 2 )]/^ - b2 


log, 


l(l - r) 


(l - r 2 )^ - b 2 r - b 2 + /(r 2 - b 2 ) (l - b 2 ) 


_ 8 _ 

3* 


(l - r2)(l - X2)3/ 2 

1(1 ~ r) 


loge 1 + 


2(1 - r 2 ) (l - b 2 ) 


3/2 


l°g e 


r - b c 


2b(r 2 - b 2 ) (l - b 2 ) 
b(l + r) 


1 - \/l - b2 


Kl ~ r) 


2r(l - r 2 ) l/r 2 - b 2 


log. 


2(l - r2)( r 2 - b2) 3 / 2 

(1 - r)(r + Vr 2 - b 2 ) 


- b 2 + \/(r 2 - b 2 ) (l - b 2 ) 


(Equation continued on next page) 
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H 3 (r,b) 


1 , -1 , (r - b 2 )(l + r) 

sin x b + * — 


2r(l - r) 


br 


2b (l - r 2 )(r 2 - b 2 )\/l - b 2 
b(l + r) 


(l - r 2 ) (r 2 - t 2 ) 3/2 
b(l - r) 


, b 
l0 Se Z + 


2(1 - r 2 ) (l - b 2 ) 


3/2 


loge 


2r(l - r 2 )/l - b2 
b(l - r) 


2r(l - i^Vr 2 - b 2 


log. 


r + 


r - b 2 + l/(r 2 - b 2 ) (l - b 2 ) 


b 1 . _1 

— + — : r sin x 


2r(l - r) 


r 2 - b 2 


2b 


4 

it 


(l - r 2 ) Vr 2 - b 2 2b (l - r 2 ) (l - b 2 ) - b' 


,2 ,4,2 _ -k2 


r~ 

b(l - 

r) £ 

4 + r 2 ) 

61 

[l - ra) 


'2 - b 2 j 

,3/2 


b(l - r) 2 (3 + r 2 ) 


3(l - r 2 )Vr 2 - b^J 


loge 


(1 - r)(r + \/r 2 - b 2 ) 


r - b 2 + ]J(t 2 - b 2 ) (l - b 2 ) 


— 

^(1 - r ) (3 

+ 

A) , 



1 

ro 

ro 

"h"' 

- b 

,)3/2 


b(l 

- r) 2 (3 + 

d 


b 

3(l 

- r 2 ) 3 /l - 

b2 

l°gg 

1 

- \/l - b2 

r 

4b(l - r) 2 



b3(l - r) 2 


3(1 - r 2 ) A - b 2 3(1 - r 2 ) (l - b 2 ) 


372 


logg b - 


(Equation continued on next page ) 


H lb 
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r - b 2 

(l " 

r ) { 2(1 - b 2 ) + (2 

- » a ) . 

6b (r^ - b 2 ) (l - b 2 ) 


6b (l - r 2 ) 2 Jl - l 

> 2 

r(l - r) (l - b 


2b(l - r) £ 

! (r ♦ b 2 ) 

6b (l - r 2 ) 2 (r2 - b 2 ' 

)/l - b : 

2 3(l - r 2 ) 2 (r 2 

- b 2 ) (l - b 2 ) 

4br(l - r) 2 


b 3 (l - r) 2 

1 np ^ + 

3(1 - r 2 ) Vr 2 - b 2 

3r(l 

- r 2 ) 2 (r 2 - b 2 ) 3/ 

/ 2 ±og e r + 

b(l - r ) 2 (l + r 

2 ) , 



6r(l - r 2 ) 2 (l - b 2 

^3/2 



b(l - r) 2 (l + 3r 2 ) 


b(l - r) 

4- 

3r(l - r 2 )Vl - b 2 

-L<Jg e 

r 

- b 2 + )/(r 2 - b 2 )(l - b 2 ) 

b(3r 2 - 2b 2 ) 


r + \/r 2 - b 2 


6r(l - r)(r2 _ b2) 3 

72 log e 

b 


(1 - r)^!- 2 - b 2 ) 

+ (r - 

b2 2 (1 - 

r)(»* - * k ) 

6b (l - r 2 ) 2 /^ 2 

- b 2 

6br(l - r 2 j 

a r 

(l - b^)\/r^ - 


/ 
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Figure 1 .- Plan form and profiles of delta wing with, varying thickness ratio. 


MCA TU 2858 






) 

i 



(a) Subsonic leading edge (b) Supersonic leading edge 

(P« < 1). (0€ > 1). 

Figure 2.- Limits of integration for deter mining the potential due to a 
source distribution representing a subsonic- and Bupersonic-leading- 
edge ving. 
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Figure 3 .- Spanwise variation of thickness distribution for various 

values of m. 
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(a) Subsonic leading edgej b = 0.2, r = 0.5, and m = - 

Figure k.- Variation of spanwise pressure distribution in tbe chordwise 
direction for single -wedge -profile, variable-thickness-ratio delta 
wings with subsonic and supersonic leading edges. 
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(c) F 3 (r,b). 

Figure 6.- Concluded 
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Figure 7.- Concluded. 
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(a) m for G- and H-f unctions . 



ft) m for F-f unctions . 

Figure 8.- Variation of m for minimum Cd/Cd ' "based on projected- 
frontal-area considerations -with 1 - r for different values of "b 
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(b) Cd/Ce' ' for F -functions. 

Figure 9-- Variation of Ce/Ce' with 1 - r for different values of b 
using the values of m from figure 8. 







Variation of Cd/Cd 1 with m for a caBe when m for 
Cd/Cd ' prescribes an impractical maximum thickness 
iion. b = 0.5, r => 0.9. 
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(a) m for G- and H-functions. 



(b) m for F -functions. 

Figure 11.- Variation of m for minimum Cd/Cd ' based on internal- 
volume considerations -with’ 1 - r for different values of b. 
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Figure 13.- Variation of Cd/Cd' with m for a case when the m for 

minimum Cd/Cd' prescribes a fictitious maximum thickness distribution. ^ 

b = 0,5, r = 0.55. 
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